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Abstract 



We apply the statistical measure of complexity, introduced by Lopez-Ruiz, Mancini and Calbet to a hard- 
sphere dilute Fermi gas whose particles interact via a repulsive hard-core potential. We employ the momentum 
distribution of this system to calculate the information entropy, the disequilibrium and the statistical com- 
plexity. We examine possible connections between the particle correlations and energy of the system with 
those information and complexity measures. The hard-sphere model serves as a test bed for concepts about 
complexity. 
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O 1 Introduction 



There has been a remarkable growth in research focused on complexity and in general on information theories in 
recent years in a variety of fields pQ including physics [3] , chemistry [3J , biology [J , neuroscience [3] , mathematics 
[6] and computer science [7]. In particular there are various applications in quantum many-body systems e.g. 
53 atoms [H [HI HO]) nuclei jTTJ (T2] [T31 [T4], atomic clusters bosonic traps [T3] e.t.c. In fact, they lead to the 
clarification of basic quantum concepts and provide results about the information content of systems according to 
various definitions e.g. Shannon information entropy [TS], Fisher [TB] information, Onicescu information energy 
[T7] e.t.c. Also they represent a suitable framework, via a probabilistic description, to assess the presence of 
interactions, correlations with experimentally measurable quantities, the derivation of universal relations e.t.c 
[n[TJ[]31[T3[T![T7l[^ Thus, traditional methods can be 

extended by an alternative information-theoretic way and give new insights for the treatment of simple quantum 
systems as well as more complicated many-body ones. A recent and important advance is to calculate several 
complexity measures, based on a probabilistic description via previous experience on information entropy in order 
to quantify statistical indicators of complex behavior in different systems scattered in a broad spectrum of fields 
[S 112 |Ml|3g 113 [3710^ Related research 

started connecting the above measures with experiment e.g. Fisher information entropy has been found to correlate 
with the ionization potential and dipole polarizability in atoms [S] and also complexity in a correlated Fermi gas 
has been connected with the specific heat Cy e.t.c. |56j . 

The statistical measure of complexity Clmc introduced by Lopez-Ruiz, Calbet and Mancini (LMC) [33J is 
defined in the form of the product Clmc = SD. Here, S is the information entropy i.e. the information content of 
the system, while D is the disequilibrium i.e. the distance to the equilibrium probability distribution. Although 
complexity is a multi-faced quantity and several definitions of complexity measures have been proposed, the LMC 
measure has been employed recently in various studies for the following reasons: it exhibits the correct asymptotic 
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properties of a well-behaved measure of complexity, as expected by intuition e.g. it vanishes for the two extreme 
cases of a perfect crystal (complete order) and ideal gas (complete disorder) . In addition it is easily calculable for a 
quantum system, described by its very nature by probability densities leading to a feasible calculation of its basic 
factors S and D. Other definitions of complexity, although sometimes may be considered that they describe certain 
aspects of complex systems in a satisfactory way, they have other disadvantages e.g. Kolmogorov's algorithmic 
complexity is hard to compute. It is defined as the length of the shortest (optimum) program needed to describe 
the system, a goal difficult to attain and prove [58] . Our approach is a pragmatic one: we start from the LMC 
definition, which is relatively easily calculable and hope to improve in the future, by a assessing the obtained 
results, comparing with other definitions of complexity e.g. the SDL measure according to Shiner, Davison and 
Landsber [37] . 

The initial definition of Clmc has been slightly modified in a suitable way by Catalan et al [34] . leading to 
the form C = e s D applicable to systems described by either discrete or continuous probability distributions. In 
[34] it was shown that the results in both, discrete and continuous cases, are consistent: extreme values of C 
are observed for distributions characterized by a peak superimposed onto a uniform sea. Moreover, C should be 
minimal, when the system reaches equipartition and the minimum value of C is attained for rectangular (uniform) 
density probabilities giving the value C = 1. Additionally, C is not an upper bounded function and can become 
infinitely large. 

The motivation of the present work, is to extend our previous study on complexity measures of uniform Fermi 
systems |13j . by employing the complexity measure proposed by Lopez-Ruiz et al. |33j . using probability dis- 
tributions in momentum space. In uniform systems the density p = N/V is a constant and the interaction of 
the particles is reflected to the momentum distribution which deviates from the theta function form of the ideal 
Fermi-gas model. Our aim is to connect C, a measure based on a probabilistic description and the shape of the 
corresponding momentum distributions to the phenomenological parameters introducing the inter-particle correla- 
tions. 

The study of uniform quantum systems (both fermionic or bosonic) in momentum space is very important. 
Very interesting phenomena such as superfluidity, superconductivity, Bose-Einstein condensation e.t.c are observed 
and also well defined in momentum space. Thus, it is interesting to concentrate our study on the connection 
between complexity, defined in momentum space to the correlated behavior of a fermionic (or bosonic ) system, by 
employing the simple, but effective, hard-sphere model. Our specific application is nuclear matter. The basic model 
with hard-spheres is a suitable starting point in order to assess the relevance of various concepts and definitions of 
complexity. The present application in a correlated Fermi system like nuclear matter is facilitated by our previous 
experience. We use the simplest potential of a hard-sphere interaction with a hard core radius. The outline of the 
present work is the following: In Sec. 2 we present the momentum distribution and information measures employed 
to quantify complexity, while Sec. 3 contains our numerical results and discussion. Finally, in Sec. 4 we exhibit our 
conclusions. 

2 Momentum distribution and information measures 

We adopt the formalism employed in our previous work [131 156] . adjusted here to a hard-sphere gas and specializing 
in nuclear matter. 

2.1 Momentum distribution 

The momentum distribution (MD) of an interacting Fermi system is given in general by the relation 



where 14 





(2) 



The normalization of n(k) obeys the relation 




(3) 



The simplest form for n(k) appears for an ideal Fermi gas. In this case n(k) is just a step function 

n Q (k) = ^-6(k F -k). 



(4) 
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The potential of the hard-sphere interaction is defined as follows, 



V(r) 



oo, for r < c 
0, for r > c, 



(5) 



where c denotes the hard core radius. 

The momentum distribution of a hard-sphere dilute Fermi gas had previously been calculated by Czyz and 
Gottfried [55] and also by Sartor and Mahaux [SD]- The above authors have studied a low density Fermi gas, whose 
particles interact via a repulsive hard core potential of the form ([5]). In this model, the quantities of interest can 
be expanded in powers of the parameter (kpc). 

The analytical expressions for the dimensionless n(k) in the hard-sphere Fermi gas model have the following 
form for k < kp [BP] 
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For y/2 < x < 3: 
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Another characteristic quantity, used as a measure of the strength of correlations of uniform Fermi systems, is 
the discontinuity, Z, of the momentum distribution at fc/fcp = 1. It is defined as 



Z = n{\-) -n(l+). 



(10) 



The behavior of the momentum distribution, as a function of a; = k/kp for various values of the correlation 
parameter kpc is shown in Fig. 1. The discontinuity Z is also displayed in each case. For ideal Fermi systems Z = 1, 
while for interacting ones Z < 1. In the limit of very strong interaction where Z — ¥ 0, there is no discontinuity in 
the momentum distribution of the system. The quantity (1 — Z) measures the ability of correlations to deplete the 
Fermi sea by exciting particles from states below it (hole states) to states above it (particle states) [61] . 

The asymptotic behavior of n(l~) for x— > 1~ reads [60j 



n(l~) w 1 
while for x— > 1 + 



3tt' 



:(v-l)(k F c) 2 



31n2 + l-3(x - l)ln | x- 1 | +(61n2- y)0- 1) 



(11) 



i(l+) « — -(y- \)(k F c) 2 31n2- 1 - 3(x - l)ln(a; - 1) - (6 In 2 - 7)(x - 1) 



(12) 



It is worthwhile to note the existence of a logarithmic singularity in the function n(k) at k = kp, a general feature 
of normal Fermi systems. 
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The discontinuity Z, according to Eqs ([TO]) and (fTT|) . (|12|) is given by 

Z = 1 - 4rln2(z/- l)(fc F c) 2 . (13) 

The energy per particle E of the ground state of v = 4-component fermion fluid of hard-spheres, in the low-density 
expansion, has been derived in Refs. [62]. Accordingly, the energy E, in units of the ideal gas energy Eq, is given 
by [62] 

E/E =e(y)~l + D 1 y + D 2 y 2 +D 3 y 3 +D 4 yHny, y = k F c : E = \%^. (14) 

5 2m 

where the coefficients Di are given in Table VI of Ref. [62] . It is one of the aims of the present work to investigate 
the connection between the various information measures and complexity with experimental quantities (as the 
ground state energy and the discontinuity Z). In addition, we intend to produce not only qualitative but mostly 
quantitative results, by connecting the strength of the correlations with the above measures. 



2.2 Information measures 

The information entropy in momentum space is given by the relation 

S k = - J n(k)\nn(k)k. (15) 
So, for an ideal Fermi gas, using Eq. (0}, Sk becomes 

S k = S = lnV k = ln(— -IV (16) 



v r. 



For correlated Fermi systems, Sk, can be found from Eq. (1151) by employing Eq. ([T]). Sk is written now 



4tt / rK 



S k =lnV k -y- k I n < {k)\n.n < {k)dk + J k 2 n > (k)lnn > {k)dk J . (17) 

The correlated entropy S k has the form 

Sk = So + S CO i, (18) 

where Sq is the uncorrelated entropy given by Eq. (|16p and .S^or is the contribution of the particle correlations to 
the entropy. That contribution can be found from the expression 

Scot = — 3 x 2 n < (x) In n < (x)dx + J x 2 n > (x) hin > (x)dx^J , (19) 

where x — k/kp. 

The disequilibrium D k (or information energy, defined by Onicescu [171). in momentum space as another 
functional of a probability distribution, in our case n(k) is given by the relation 

D k = J n 2 (k)dk. (20) 

For an ideal Fermi gas, using Eq. (j4]), becomes 

D k =D = ^-. (21) 
V k 

In the case of correlated Fermi systems, D k is written as 



1 4tt 



D k = TFTr\ I k 2 n 2 < (k)dk+ k 2 nl(k)dk). (22) 

Vk Vk yJo J k + J 

The correlated disequilibrium D k is 

D k = D Q D COI , (23) 
where Dq is given in Eq. (|2ip and D cor can be found also from the expression 

D cor = iij x 2 n 2 < (x)dx + x 2 n 2 > (x)dx\ . (24) 
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The statistical complexity measure, proposed by Catalan et. al. (34] . in momentum space, is denned as 

Clmc = C = DkHk, (25) 
where H represents the information content of the system defined as 

H k =e s «. (26) 

It is easy to show that 

C = C a C cm = D cor e s ™, C = D e s ° = 1. (27) 

The physical meaning of Eq. (|2"T)) is clear. In the case of an ideal Fermi gas (see Eq. (0|) C is minimal with 
the value Co = 1 (see also [34]). Moreover as pointed out in Ref. [34], C is not an upper bounded function and 
can therefore become infinitely large. From the above analysis it is clear that complexity C is an accounter of 
correlations in an infinite Fermi system. Thus, the next step is to try to find the connection between C and the 
correlation parameters of the systems. The correlation invoke diffusion of the momentum distribution and we 
expect this effect to be reflected on the values of C. 



3 Results and discussion 

The behavior of the momentum distribution, as a function of fc/fc^, for various values of the wound parameter fej?c 
is shown in Fig. 1. The discontinuity Z is also displayed in each case. For ideal Fermi systems Z = 1, while for 
interacting ones Z < 1. In the limit of very strong interaction Z — 0, there is no discontinuity on the momentum 
distribution of the system. 

The calculated values of 5 cor , D cor and C for nuclear matter versus the correlation parameter kpc are displayed 
in Fig. 2. S COI and C increase with y, while D cor decreases. We fitted the numerical values of the above quantities, 
with simple functions of y = UfC and we find respectively the following formulae 

S cor (y) = ay , a = 2.16379, /3 = 1.67053. (28) 

D cor (y) = 1 + a/, a = -0.79871, /3 = 1.83155. (29) 

C{y) = 1 + a/, a = 1.68358, /3 = 1.67566. (30) 

The values of the parameters a, (3 and 7, for each case, have been selected by a least squares fit (LSF) method. It 
is worthwhile to notice that for an ideal Fermi gas there is an upper limit for C (C max ~ 1.5845) [55] . However, for 
an interacting Fermi gas there is no such constraint (at least to the region under consideration in the present work). 
It is obvious that the calculated values of complexity reflect the different way that the interaction or temperature 
affect the trend of the momentum distribution. 

The quantity (1 — Z) measures the ability of correlations to deplete the Fermi sea by exciting particles from 
states below it (hole states) to states above it (particle states) [51] , The dependence of S cm , D cm and C on the 
quantity (1 — Z) is shown in Fig. 3. It is seen that S cm and C are increasing functions of (1 — Z), while D COI is 
a decreasing one, as a direct consequence of the dependence of the above quantities on the correlation parameter 
kpc. That dependence can be reproduced very well by simple expressions as in Eqs. (j2"5|) . ([2U]) and (IBTfl) replacing 
fcpc by (1 — Z) 

S cor (Z) = a(l - Zf, a = 2.49614, j3 = 0.83527. (31) 
D cm {Z) = 1 + a(l - Zf, a = -0.93413, /3 = 0.91574. (32) 
C(Z) = 1 + a(l - Z) fi , a = 1.94305, /3 = 0.83784. (33) 

From the above analysis we can conclude that LMC complexity C can be employed as a measure of the strength 
of correlations in the same way the wound and the discontinuity parameters are used. An explanation of the above 
behavior of C is the following: The effect of nucleon correlations is the departure from the step function form of 
the momentum distribution (ideal Fermi gas) to the one with long tail behavior for k > hp. The diffusion of the 
distribution leads to a decrease of the order of the system (the disequilibrium Z\. decreases and the information 
entropy Sk increases accordingly). In total, the contribution of Sk in C dominates over the contribution of Dk and 
thus the complexity increases with the correlations (at least in the region under consideration). 

It is one of the aims of the present work to connect C, a measure of complexity based on a probabilistic 
description and the shape of the corresponding momentum distributions with other quantities like the ground state 
energy per particle E. In view of the above, we display in Fig. 3 the dependence of S cor , D cor and C, as well as the 



5 



energy fraction e(y), on the correlation parameter Ufc. The dependence of S CO r, D cor and C on e(y), as displayed 
in Fig. 4, is in a very good approximation linear. The fitted expressions are the following: 

e = 1.0479 + 1.273535^, (34) 

e = 4.6861 - 3.5985L> cor , (35) 

c = -0.5847+ 1.6358C. (36) 

In total we observe an empirical connection of the energy with S cor , D cor and C calculated employing information 
entropy which, by definition, is not related directly to the energy of the system, in contrast to the traditional 
concept of thermodynamic entropy. The above results confirm our recent finding, according to similar lines, that 
there is also a connection between the "energy-like" quantity specific heat Cy of an ideal electron gas with the 
complexity C [5B] . 

Finally, in Fig. 5 we compare the present results with those taken from the Low Order Approximation (LOA) 
method. Thus the momentum distribution of nuclear matter is evaluated by employing the LOA and the MD takes 
the form |61j 

n h0 A(k)=9(k F -k)[l-k dir + Y(k,8)] + 8[k dil Y(k,2)-[Y(k,A)} 2 ], (37) 



Q-k+- e -k- r k + , . r\ k -\ 



2k 



c- x Y{k, n) = -= + / e~y y + sgn(fc_) / e"» y, (38) 



o Jo 



1 //i\ 3 / 2 r k ~ kp ± k , . 

= — =' L i • k=——, k± = , n = 2,4,8. (39) 



where 



and 

while sgn(x) — xj \ x |. The dimensionless wound parameter fcdir can serve as a rough measure of correlations and 
the rate of convergence of the cluster expansion is defined as 

k dil =pjlf(r)-l] 2 r. (40) 

The normalization condition for the momentum distribution is 

n LOA (fc)fc 2 k = - k F . (41) 

o 

The following relation between the wound parameter fcdir and the correlation parameter /3 

1 Af x3 



^'urAf)- (42) 

It is clear that large values of fcdir imply strong correlations and simultaneously poor convergence of the cluster 
expansion. In the numerical calculations the correlation parameter f3 is in the interval: 1.01 < (3 < 2.482. That 
range corresponds to 0.3 > fcdir > 0.02 and this is reasonable, in the case of nuclear matter [61]. However, the origin 
of the two methods is different and as a consequence they influence in a different way the momentum distribution. 
Nevertheless, we found that S cor , D cor and C exhibit a similar trend as a functions of the parameter 1 — Z . 



4 Conlusions 

In conclusion, we calculate information and complexity measures of a uniform Fermi system, like nuclear matter, 
in the framework of the hard-sphere model. The effect of correlations is connected intuitively with the concept 
of complexity, in a qualitative, and somehow vague way as stated in [56j as well. In fact, it turns out that all 
information measures used by us, show a strong dependence on the correlation parameter kpc as well as on the 
Fermi discontinuity Z. The most distinctive feature is the occurrence of an almost linear dependence between 
information measures and complexity on energy. The above statement is in keeping with the recent finding of the 
existence of an empirical connection between the specific heat and complexity [56] . However, the applicability of 
our approach is much wider than nuclear matter, since the impenetrable hard spheres (not overlapping in space) 
can simulate the extremely strong repulsion that atoms and spherical molecules feel at very small distances. Thus, 
the significance of a suitable quantification of complexity emerges in statistical mechanics of fluids and solids. We 
do not claim that our approach is the only or more important one, but so far our results are interesting and 
encouraging. We stress again that the proposed LMC measure of complexity is by its definition an appropriate one 
and specifically tailored for systems described statistically, through a probability distribution. Since "information 
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is physical" [63], it is promising to examine how far this quotation goes, in the sense that calculations employing 
a good measure of information content of a quantum system (and consequently of complexity) are expected to 
give, at least, interesting results of physical relevance. Landauer dedicated his research on similar ideas. Hence, 
in our present work and previous ones, we proceed towards a numerical quantification of complexity. One of our 
goals is to examine, as a first step, whether a particular definition of complexity is reasonable and robust enough 
e.g. if one increases the value of one parameter (or parameters) describing a quantum system, the corresponding 
value of complexity increases accordingly. This would satisfy minimally intuition about complexity and validate 
its definition. 
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Figure 1: (a) The momentum distribution for correlated nuclear matter versus k/kF for various values of the 
correlation parameter kpc (b) The same, but on a logarithmic scale, for n(x). 
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Figure 2: (a) Scot, D cor and C of nuclear matter versus the the correlation parameter kpc and (b) the discontinuity 
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Figure 3: S cor , D cor , C and the energy fraction e (for definition sec Eq. (fT^f) versus the correlation parameter fcpc. 
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Figure 5: S CO r, D cor and C of nuclear matter versus the the discontinuity parameter (1 — Z) calculated by employing 
two different methods: the Low Order Approximation (LOA) (see text) and the hard sphere interaction (HS). 
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